The inverse cascade of magnetic helicity in magnetohydrodynamic turbulence 
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The nonlinear dynamics of magnetic helicity, H , which is responsible for large-scale magnetic 
structure formation in electrically conducting turbulent media is investigated in forced and decaying 
three-dimensional magnetohydrodynamic turbulence. This is done with the help of high resolution 
direct numerical simulations and statistical closure theory. The numerically observed spectral scaling 
of H M is at variance with earlier work using a statistical closure model [Pouquet et al., J. Fluid 
Mech. 77 321 (1976)]. By revisiting this theory a universal dynamical balance relation is found 
that includes effects of kinetic helicity, as well as kinetic and magnetic energy on the inverse cascade 
of H M and explains the above-mentioned discrepancy. Considering the result in the context of 
mean-field dynamo theory suggests a nonlinear modification of the a-dynamo effect important in 
the context of magnetic field excitation in turbulent plasmas. 
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The emergence of large-scale magnetic structures in 
turbulent plasmas is dynamically important in many as- 
trophysical settings, e.g. with regard to the interstellar 
medium or the magnetic field generation in planets and 
stars by the turbulent dynamo effect, e.g. [3-Q- The 
structure formation can be studied via the magnetic he- 
licity H M = w Jv dVa ■ b where b is the magnetic field 
and a denotes the magnetic vector potential. This topo- 
logical characteristic of magnetic fields yields a measure 
of the linkage and the twist of the field lines, @, @. In 
the magnetohydrodynamic (MHD) single-fluid approxi- 
mation [fjj which neglects microscopic scales and the as- 
sociated kinetic dynamics H M is ideally conserved in a 
three-dimensional volume with periodic or closed bound- 
ary conditions Q- It is thus prone to a nonlinear and 
conservative inverse spectral cascade process in the iner- 
tial range of MHD plasma turbulence. If driven at small 
scales, £, the cascade results in spectral transfer of mag- 
netic helicity towards small spatial wavenumbers k ~ 
Q, i.e. to the formation of large-scale magnetic struc- 
tures. This process is thus of fundamental importance 
with regard to, e.g., the dynamics of magnetic fields in 
the above-mentioned turbulent astrophysical settings. In 
spite of its importance little is known about the nonlin- 
ear dynamics which underlies the inverse cascade. It is 
the purpose of this work to shed some light on the rather 
mysterious nonlinear phenomenon creating large-scale or- 
der out of quasi-random turbulent magnetic fluctuations. 
Please note that the constraining effect of magnetic helic- 
ity conservation on certain a-dynamo configurations is, 
although important, beyond the scope of this work, cf., 
for example, |9j]. The existence of such an inverse cascade 
was first demonstrated in numerical simulations based 
on the eddy damped quasi normal Markovian (EDQNM) 
closure model of three-dimensional MHD turbulence [10| 
which is to the best of our knowledge the only work 
dealing theoretically with the spectral self-similarity of 
magnetic helicity. The associated self-similar spectral 
signature in the turbulent inertial range, ~ fc~ 2 , is in 



agreement with dimensional analysis based on a con- 
stant nonlinear spectral flux. Several studies applying 
direct numerical simulations (DNS) find inverse transfer 
of magnetic helicity, see e.g. [Ill4l3j . without reporting 
self-similar scaling behaviour, a notable exception being 
13 [ . 

In this Rapid Communication the inverse cascade of 
magnetic helicity in homogeneous MHD turbulence is 
studied by three-dimensional high-resolution direct nu- 
merical simulations. In the main setup kinetic and mag- 
netic energy and magnetic helicity are injected at small 
scales of the initially excited spectral range of turbulent 
fluctuations. It is shown that macroscopic quantities, 
in particular kinetic helicity and the ratio of kinetic and 
magnetic energy, have an important influence on the H M - 
cascade that is captured by a universal relation based 
on dimensional analysis of the MHD-EDQNM closure 
model. This insight suggests a new link between mag- 
netic helicity and mean-field dynamo theory in particular 
with regard to the saturation behaviour of the dynamo 
mechanism. 

The dimensionless MHD equations are written as 

d t u = Vx(vxw+jxb) + (^ n (-l)"- 1 A n +AA- 2 )a;+F v , 

(1) 

ftb = Vx(vxb) + (77n(-l) n_1 A" + AA- 2 )b + F b , (2) 



V b = 0, 



(3) 



where v is the velocity, u> = V x v the vorticity, b the 
magnetic field, and j = V x b the electric current density. 
Equations CO)-® are solved by a standard pseudospec- 
tral method using a leapfrog scheme for time integra- 
tion. Anti-aliasing is achieved by spherical mode trun- 
cation. The simulation domain is a triply 27r-periodic 
cube discretized by 1024 3 collocation points. Hypervis- 
cous small-scale dissipation operators of order n — 8 are 
used to improve scale separation parametrized by the hy- 
pcrdiffusion coefficients \x n and rj n with the hyperviscos- 
ity fig = 9 • 10 -41 and ^s/vs — 1- Boundary effects at 



FIG. 1: Inverse cascade of mag netic helicity, Hjr, in 3D-MHD 
turbulence for t 6 [1, 6.66]. Curves are given for points in time 
spread equidistantly over the interval. Magnetic fluctuations 
with maximal magnetic helicity and non-helical velocity fluc- 
tuations are supplied through the forcing with fco = 206. 

smallest spatial wavenumbers fc are alleviated by a large- 
scale energy sink AA -2 for both fields with the constant 
A set to 0.5. The forcing terms F v and Fb are ran- 
dom, delta-correlated processes of equal amplitude that 
act over a band of wave numbers fc 6 [kg — 3, fco + 3] with 
fco = 206. They inject velocity- and magnetic-field fluc- 
tuations with well defined magnetic and kinetic helicity, 
kinetic helicity being defined as H K — ^7 J v dVv ■ u>. 
Such driving, chosen here for simplicity and numerical 
efficiency, could in principle be realized by a random 
small-scale distribution of electric currents and forces. 
The initial velocity and magnetic fields are smooth with 
equal energies, random phases and fluctuations that have 
a Gaussian energy distribution, peaked at fco- In the 
course of the simulation the total energy quickly attains 
a quasi-stationary state, fluctuating around unity with 
E M /E K « 9. Cross-helicity, H c = ±, J y dVv ■ b, is 
negligible. The simulation is carried up to t=6.66 large- 
eddy turnover times. The application of hyperviscous 
dissipation operators while necessary to observe well- 
developed scaling ranges precludes the unambiguous def- 
inition of a Reynolds number. 

The temporal evolution of the magnetic helicity spec- 
trum over the simulation period shown in Fig. [I] indi- 
cates inverse spectral transfer. This is also reflected by 
the nonlinear spectral flux of magnetic helicity, II? = 
2 J^dk'k' 2 J dtt(bl, ■ [v x b] k < +c.c.) with '[.] k ' denot- 
ing Fourier transformation and '*' standing for com- 
plex conjugate (c.c). The flux spectrum shown in Fig. 
[2]Ja), is constant over finite wavenumber intervals on 
both sides of the forcing band signalling equilibrium 
between turbulence driving and dissipation. While on 
the right-hand-side direct spectral transfer is observed 
as a result of the small-scale energy sink, on the left- 
hand-side of the forcing band an inverse cascade de- 
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FIG. 2: Moduli of fluxes of magnetic helicity (a) and total 
energy (b) (inverse: solid curves, direct: dotted curves) nor- 
malized by corresponding dissipation rates. Flux spectra are 
shown for t 6 [1.66, 6.66] for the forced simulation and are 
spread equidistantly over this period. 

velops which is driven by the constant magnetic helic- 
ity input around fco- The spectral flux of total energy, 
E = E K + E M = Jd 3 fc(|w k | 2 + |6k| 2 )/2, is given by 
n fc = J fc dk'k 12 J dQ(ik' x [v x w + j x b] k , ■ v*, + ik' x 
[v x b] k < • b k , + c.c.) and is shown in Fig. HJb). It lacks 
spectral constancy in the inverse cascade region and it 
is principally carried by magnetic energy transfer. This 
suggests that the inverse energy flux is a consequence 
of the inverse cascade of H M . The linear scaling of the 
life-envelope which follows from the dimensional estimate 
Ek ~ kH^f in combination with the approximate con- 
stant value of the envelope of 11^ also supports this 
interpretation. 

The compensated magnetic helicity spectrum at the 
end of the simulation period is displayed in Fig. [3ja). 
It exhibits two approximate scaling ranges: on the di- 
rect transfer side for 250 < k < 400 and in the inverse 
transfer region, 7 < fc < 30. The corresponding asymp- 
totic scaling laws are H¥- ~ fc -3 3 (inverse, cf. [13j]) and 
H^f ~ fc~ 15 (direct). The latter value has to be taken 
with care due to the very short spectral range and the 
high-order hyperviscosity acting at largest wavenumbers. 
The inverse cascade scaling is at variance with the fc -2 - 
behaviour reported in (loj . 

bmcc H M is the helicity of the magnetic vector poten- 
tial, a, its inverse cascade pulls quantities derived from 
this potential, e.g. magnetic energy and to a lesser extent 
also electric current density, jfc ~ fc 2 afe, towards large 
scales. Magnetic and velocity field are intrinsically cou- 
pled in MHD turbulence by Alfvenic fluctuations [14[ and 
thus similar behaviour is observed for the spectral kinetic 
energy and kinetic helicity, H^, as well. The respective 
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FIG. 3: Compensated spectra of (a) magnetic helicity , (b) to- 
tal (solid line), kinetic (dashed line), magnetic (dashed-dotted 
line) energy, (c) kinetic helicity. (d) A = {E* / E™y< H J k / 
for 7=1 (dashed line) and 7 = 2 at t=6.66 in the driven 
simulation. 



i.e. E k ~ k~ 5 ' 3 and Eg = \Ejf - E£\ ~ k' 1 1 3 (not 
shown) with an excess of magnetic energy. The obser- 
vations are in agreement with the interpretation of the 
finite levels of E^f and of H k as a result of the local and 
temporary stirring induced by changes of magnetic-field 
topology. This is to be expected in the course of the in- 
verse cascade of magnetic helicity. As will be shown in 
the following, the lacking equipartition of Ef? and E^f 
and the presence of kinetic helicity are the reasons of 
disagreement with the above-mentioned EDQNM simu- 
lations of fioj . There, the relaxation time of nonlinear in- 
teraction, Okpq, which represents a free parameter of the 
EDQNM approach and determines the nonlinear process 
governing turbulent dynamics, is chosen to be the Alfven 
time, (fci?o) _1 - Consequently, the resulting dominance 
of Alfvenic interactions drives the system in the iner- 
tial range into nearly perfect equipartition of kinetic and 
magnetic energies. 

The observed spectral scaling of magnetic helicity 
can be better understood with the help of the integro- 
differential EDQNM equation for the evolution of H^f . A 
formally similar approach has been successful with regard 



to the residual energy spectrum Ej? = \Ejf — E*\ [14 1. 
The equation for the evolution of Hfc 1 in the EDQNM 
model [10( reads 



spectra are thus expected to inherit the self-similarity 
from H]f as is indeed observed, see, e.g., Fig. EJc). It 
seems to be reasonable to regard the merging of current- 
carrying structures where the currents have significant 
positively aligned components and thus experience mu- 
tual attraction as the physical mechanism of the inverse 
cascade of magnetic helicity, cf. [15j . This is also in 
line with recent statements in the literature about non- 
locality of the magnetic helicity flux 12, 3]. Terming the 



inverse transfer of H^r a "cascade" is thus merely follow- 
ing convention and not a physical description of the ac- 
tual process of a spectrally non-local merging of current- 
carrying strucures. Note that in this simulation no ki- 
netic helicity is injected by the turbulence driving. In 
the case of simultaneous injection of kinetic helicity the 
spectral diagnostics discussed in this work show no sig- 
nificant difference to the presented simulations. In gen- 
eral the details of the small-scale forcing, like random- 
ness of amplitude and/or complex phases, were verified 
to have no measureable impact on the self-similar be- 
havior reported below. The only significant parameter 
in this respect is the separation of the smallest admis- 
sible wavenumber, k = 1, and the forcing- wavenumber 
of about two decades. The scale-separation determines 
the extent of the observable self-similar inverse-cascade 
dynamics. It is thus even more important than the classi- 
cal Reynolds number which is not well defined due to the 
necessary high-order hyperviscous small-scale diffusion. 

The energy spectra in Fig. EJb) exhibit approximate 
scaling known from decaying turbulence, cf., e.g., 14j.ll7j|. 



{d t +r n k i )H?= dpdq9 kpq (Tt v + Tr + Tn (4) 
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1 k - e kpq—E k H q -]kpq—E q H k . 

The geometric coefficients h kpq , e k p q , and j kpq follow 
from the solenoidality constraints Q and are given in 
[Tol | ■ T ne 'A' restricts integration to wave vectors k, p, 
q which form a triangle, i.e. to a domain in the p-q plane 
which is defined by q = p + k|. The time 9 kpq is char- 
acteristic of the eddy damping of the nonlinear energy 
flux involving wave numbers k, p, and q. It is defined 
phenomenologically but its particular form does not play 
a role in the following arguments. 

The three nonlinear contributions on the right-hand- 
side of Eq. (jH) can be associated with the advective 
(T| dv ) and explicitly twisting (T^ hl ) effects of turbulent 
fluctuations, as well as self-interaction (T^ ot ) of the mag- 
netic field through the Lorentz force. Assuming that the 
most important nonlinearities involve the turbulent ve- 
locity and that the spectral scaling range of Hjf is sta- 
tionary, a dynamical equilibrium of turbulent advection 
and the i? M -increasing effect of helical fluctuations is as- 
sumed, i.e. Tt dv ~ Tfc hl . Dimensional approximation of 
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the respective flux terms, kEg H^ 1 ~ k 1 E£ I Hg , yields 

H?~(^yh?H?, 7 = 1- (5) 

This is a statement about the spectral dynamics of kinetic 
and magnetic helicities (or kinetic and current helicities, 
H J = 2^ J V {V x b) • b, since H( ~ k 2 H^) in the case 
of Eg/E^f ^ 1. The agreement of Rel. © with the 
numerical experiment is significantly improved, see Fig. 
[3f d) , when increasing 7 by one. This yields the main 
result 

Rel. © is fulfilled, i.e. constant and close to unity, for 
almost all wavenumbers k > 12 excluding the drive and 
deep dissipation scales. It does however not cover the full 
spectral scaling range of H^ 1 due to its susceptibility to 
the asymmetry of energies and helicities introduced by 
the large-scale energy sink. The higher-order modifica- 
tion of Rel. ([5]) can not be motivated within the frame- 
work of quasi-normal EDQNM theory using the chosen 
approach of nonlinear equilibrium. The underlying cause 
is presumably the non-locality of the inverse cascade pro- 
cess which is not captured by the dimensional simplifica- 
tion of the EDQNM equations. 

For verification purposes a test simulation of decaying 
turbulence with the same numerical resolution of 1024 3 
and an initially finite level of magnetic helicity, i.e. 50% 
of the energetically possible maximum H^f ax ~ E M /ko, 
is conducted. The initial Gaussian random-phase energy 
spectrum with Eg — Ej^ 1 for all k is peaked around 
fco = 70 to allow some development of inverse trans- 
fer. The distribution of Hjf is homogeneous over the 
initial spectrum. The hyperdiffusive coefficients are cho- 
sen as = ??8 = 3 • 10~ 41 with forcing switched off. Fig. 
2] displays the most important results of this decaying 
test simulation after about 9.2 large-eddy turnover times. 
The magnetic helicity, Fig. EJa), exhibits approximate 
self-similar spectral scaling ~ k~ 3 e with most of the ex- 
cited scales in the range 15 < k < 60 having developed 
during the decay. Kinetic helicity, see Fig. |U[b), reflects 
this as it is not injected into the system but generated 
via the Lorentz force during the inverse transfer of H^f . 
It shows two distinct asymptotic power-laws: ~ k~ 5 for 
60 < k < 120 and - k 0A for 11 < k < 45 in the range 
covered by the inverse transfer of H^ 1 during turbulence 
decay. The signatures of kinetic, magnetic and total en- 
ergies are similar to the observation made in the driven 
case. The spectral flux of magnetic helicity, Fig. Etc), 
exhibits a split near ko as for k > ko small-scale dis- 
sipation determines the helicity transfer direction while 
for k < ko an inverse transfer with fc-dependent II? is 
observed. This is reflected by Rel. (JB]), cf. Fig. E£d), 
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FIG. 4: Compensated spectra of (a) magnetic helicity , (b) 
kinetic helicity. (c) Spectral magnetic helicity flux normalized 
by corresponding dissipation rate (Inverse: solid curve, direct: 
dashed curve), (d) A = (Eg / 'E^ 1 ) 1 H( / Hg for 7 = 1 and 
scaled by a factor 1/3 (dashed) and for 7 = 2. All spectra 
shown have been taken at t = 9.2 in the decaying simulation. 

which is only fulfilled in the region of approximately con- 
stant 11^ and is apparently independent of the transfer 
direction. 

Comparison of the presented findings with mean-field 
dynamo theory, cf., for example, [3j and references 
therein, in particular the a-coefficient including the cur- 
rent helicity, a ~ H J — H K , is interesting. Rel. © 
is consistent with vanishing a since the cascade process 
does not generate magnetic flux and, thus, does not act 
as a dynamo itself. It furthermore suggests a modified 
a ~ (E K / E M ) 2 H J — H K which includes the squared ra- 
tio of kinetic and magnetic energy ~ (E M )~ 2 as a turbu- 
lent dynamo-quenching mechanism of the current helicity 
contribution to a. This form of a-quenching has recently 
been observed in a numerical test-field model j 1 81 ] - 

In summary a new and probably universal relation con- 
necting the spectral behaviour of magnetic and kinetic 
helicities and energies in homogeneous MHD turbulence 
is found motivated by statistical EDQNM closure results. 
DNS of MHD turbulence that is decaying or driven at 
small scales confirm the validity of the findings for spec- 
tral intervals of constant flux of magnetic helicity. The re- 
sult has interesting connections to the a coefficient known 
from mean-field dynamo theory: it suggests an inherent 
and strong quenching of the dynamo, in particular of the 
current helicity effect, if the energy of turbulent magnetic 
fluctuations grows compared to the kinetic contribution 
of velocity. 

The authors are grateful for discussions with A. Bran- 
denburg, P. Diamond, U. Frisch, R. Grappin, P. Mininni 
and A. Pouquet. 
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